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Abstract 

The method of zero-range potentials is generalized to account for the molecular electron 
excitation process. It is made by a matrix formulation in which a state vector components 
are associated with a scattering channel. The multi-center target is considered and the 
model is applied to the example of e + Hi low energy scattering. The results of evaluation 
of cross-sections are compared with ones of the MCF and SMC methods. 



1 Introduction. 

The ideas of zero range potential (ZRP) approach were recently developed to widen limits of 
the traditional treatment Q, the author accounts higher momentum partial amplitudes || for 
a multi-center problem. Historically, the ZRP notion was introduced, perhaps, in || to model 
a potential action sustaining some parameters values. The simplification of the theory is such 
that a problem of differential equation goes down to some algebraic one M. There is the book 



and the review |L0| in which the theory is comprehensively described. 

The target of this paper relates to the other limitation of the theory connected with multi- 
channel character of real scattering phenomenon. Following the ideas of || we introduce a 
vector formulation of the scattering problem with the components that stands for different 
possibilities (channels) of the process to be described. The scattering amplitude hence is a 
matrix connecting asymptotic states. In the Sec. 1 the formalism is introduced for a multi- 
center problem. 

As an important example we consider applications to a diatomic molecule. In the framework 
of the generalization to be introduced we revisit results of the papers M, M starting from 
adiabatic nuclei problem (Sec. 2). The rotations of a molecule are included similarly to the 



mentioned approach [^] but with some difference in the formalism. Oscillations are described by 
means of the Morse potential with eigen functions - vibrational harmonics that are proportional 
to Laguerre polynomials (Appendix). The results of differential and integral cross-sections 
evaluation are compared with the direct numerical solution by some standard molecular orbits 
model (Sec.3). 



2 The matrix zero-range potentials 



As it is known, if we diminish the range and at the same time increase the depth of a pertinent 
spherically symmetrical potential then the only characteristic parameter of a potential may 
remain fixed (usually it's a scattering length or a position of the bound-state energy level), 
others characteristic are lost. This procedure allow to replace the potential action by the 
boundary condition on the wave function ^(r) at the potential center (see Ref. |], |fl) 
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where atj - inverse scattering length, Rj - position vector of the potential center j. It's necessary 
to remark that ZRP can be introduced in another ways. For example, it can be defined by the 
following equality (see Ref. ||, |§) 
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where S(pj) - Dirac function in the three-dimensional space. In order to adapt ZRP model for 
multichannel scattering we will replace the parameters atj by matrices Aj and a wave function 
ip(r) by the vector function \P(r) 
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The components ip n ( r ) of the vector function \I/(r) have the form (see Ref. [l|]) 
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(2.1) 



Here Cj are some constant vectors, k n - electron's momenta in the channel n, ko - momentum 
of the incoming electron. The representation (|2.1|) indicates the component ip n {^) describes the 
electron scattering in the channel n. The boundary conditions for components ip n { r ) have the 
form 
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This equation is basic in the ZRP theory. General properties of the matrices Aj are determined 
by symmetry of the quantum states. It means the matrices of different centers are linked by 
symmetry transformation. For example, in the case of only two centers the matrix A 2 is a 
transform of A\. 
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The diatomic, homogeneous molecules at the N-H-state approximation: Let us consider 
the problem of electron-impact excitation of a diatomic, homogeneous molecule in the iV-state 
model. The coordinate systems will have the common origin at the center of mass of a molecule. 
Then R x = — R and R2 = R, here 2R is the internuclear distance (see Fig. |IJ). The space- fixed 
axis oz of the first frame is chosen along initial momentum ko (the so-called LAB-frame, see 
the review ||10||), whereas axis oz' of the second system we direct along a symmetry axis of the 
molecule (i.e. along R, such coordinate system is known as BODY-frame). 

Though atoms in the molecule are identical the off-diagonal elements of the atom potentials 
Vi t 2 (and therefore matrices A^) can be different in sign because the parities of the molecular 
states \n) can be different. Suppose hence 

At = diag(a , . . . ,a N ) + A, (2.3) 

if A is some Hermitian matrix with zero in-diagonal elements. Let us consider the diagonal 
matrix a v = diag(r/o, • • • ,Vn), where rj n are parities of molecular states |n). Then the matrix 
A 2 is given by 

A 2 = a v A 1 a v = diag(a , . . . , a N ) + a v A a rp (2.4) 

where a n are some parameters of the molecular states \n), A - matrix of the coupling channel 
parameters. 

If the nuclei are held space- fixed then the components ip n {^) of the vector function \l/(r) 
have the form 

Mr) = e^ '^ 5 n0 + (C x ) n + (C 2 ) n , p lj2 = |r±R|. (2.5) 

Pi P2 

Therefore transition amplitudes in fixed-nuclei approximation can be constructed by the ex- 
pressions 

/„(k,k 0j R) = (CO^fc 11 ) + (C 2 ) n e- z ( k ' R ), k = h n - } (2.6) 

r 

where k - outgoing electron momentum. 

It is necessary to note that other notations are more convenient 

(<?!)„ = 5« + St\ (C 2 ) n = VoVn - 5H) . (2.7) 

Using the boundary conditions ( |2.2|) we obtain the vectors and S^~\ To introduce the 
matrices 

A( ± )=diag(^ ±) ,...,4 ±) ) + A, 
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where we used the following notation 

^2ik n R 

= a n + ik n ± r] r] n ^ R , 

one obtains for and the expressions 

= - cos(k , R) (AW-%, = i sin(k , R) (2.8) 

Using Eqs. ( |2.6|) , (|2.7|) and ( |2.8| ) we arrive to the observation that if r/o^n = +1 then the 
transition amplitude / n (k, k , R) is given by 

/ n (k, k , R) = -2 cos(k, R) cos(k , R) - 2 sin(k, R) sin(k , R), (2.9) 

otherwise if r/o^n = —1 then the amplitude / n (k, ko,R) have another form 

/ n (k, k , R) = 2z fij^ cos(k, R) sin(k , R) - 2i ^ +) sin(k, R) cos(k , R), (2.10) 



where factors = (A^ 1 ) n o- The first item of the Eqs. ( |2.9| ) and ( |2.10| ) is even function 



(under reflection k — > — k) that conform to the wave and the second item is odd function 
that corresponds wave. 

The example of diatomic, homonuclear molecules in the 2-H-state approximation: it's a 
special case of the preceding model and so we describe the case briefly, introducing additional 
useful notations. The matrices Ai and A 2 are given 

A!=( ? A V A 2 =( a \ (2.11) 

where A is the real coupling channel parameter. The equations both for amplitudes / n (k, ko, R) 
(see Eqs. (|2.9|) , (|2.10|) ) and for factors Vt^ remain valid in this case. In particular the factors 
have the form 

( ^ \ = 1 ( o{ ±] \ 
{ ni ±J ) o^eP - a 2 v -A / 



(2.12) 



3 The adiabatic-nuclei approximation 

Further we suppose that adiabatic-nuclei approximation is valid for both rotations and vibra- 



tions account. Initially this approximation was applied by Drozdov [11], Chase [12] and Oksyuk 



13 1 . The adiabatic approximation in ZRP model was developed by Demkov and Ostrovsky [|IJ 



and Drukarev and Yurova |7j (see also Ostrovsky and Ustimov [|TJJ). Differential cross sections 



of the electron-rotational-vibrational transitions can be expressed via corresponding matrix el- 
ements of the electron transition amplitude which obtained in space-fixed nuclei approximation 
by the formula 

^-(nv'j'm' «- Ovjm) = | (nv'j'm\f n (k, k , R)\0vjm ) | 2 , 
a\l k 
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where (Ovjm) - initial quantum numbers and (nv'j'm') - final quantum numbers of the electron- 
rotational-vibrational molecular states 



[R\nvjm) =R- 1 X nv (R)YJ n (R). 



Here l^ m (R) are spherical harmonics [15| and X nv (R) are vibrational harmonics (see Appendix). 
Further, we omit 'A' in arguments if it does not lead to misunderstanding. In general, the 
vibrational molecule energies essentially exceed the rotational energies thereby we can neglect 
the rotational contributions to the vibrational harmonics. The outgoing electron momentum 
k n = |k| and incoming electron momentum ko = |ko| are directly related by law of conservation 
of energy (we use atomic units throughout) 

P P 

n _L TP — _L TP 
— + &nv' — — + &0v, 

where E nv > - electron-vibrational state energy, some function of the arguments n, v' . 

Since the rotational energy levels of a molecule are so closely spaced, it is often impossible 
to resolve particular final states in a scattering experiment. In this case observed differential 
cross sections are averaged over initial rotational states (j, m) and hence it is summed over the 
final rotational states (f, m') 

da 00 b ■ j 00 j ' da 

—{nv' ^0v) = J2 tt^-t E E E -7^{nv'l'm <- Ovlm), 

Uii j=0 J m=-j j'=0 m'=-j' 

bj are relative populations of the rotational states. The summation over initial and final ro- 
tational molecular states can be exactly realized by using some simple properties of spherical 
harmonics. Final formula for observed differential cross sections have the form 



/7/j 1 U r ^ poo 2 

' W^0u) = — / X nv ,(R)f n (k,k ,R)X 0v (R)dR 

Air kn J Jo 



(3.13) 



_( n ^0 W ) = -^y ^iXo^^^ydRl/^k,^,^! 2 . (3.14) 



dVt 47T kn 

Differential cross section for pure electronic transition can be obtained by summation over 
complete set of final vibrational substates (generally - including continuous spectrum) 

da , 1 k n 

An k Jo 

Rigorously speaking, the formula is certainly correct only if the electron-vibration state energy 
E nv i and, therefore, outgoing electron momentum k n are independent on vibrational quantum 
number v ' . 

Integral cross section of the electronic-vibrational and pure vibrational transitions can be 
evaluated by integration over all angles k 

a{nv' «- Ov) = [ dk ^r(ra/ «- Ov). (3.15) 



dtt 

Differential cross sections of pure electronic transitions: Let us start from a differential 
cross section for pure electron transitions. Taking into account the analytical form Eqs. ( |2.9| ) 
and (|2.10|) we can perform exact averaging over all molecular orientations (see Eq. ( p.!4| )). 



Final result is given by 

|kn<-0i;) = £ [°°dR \X 0v (R)\ 
d\i kn Jo 
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fiW 2 „4 (+) (k, ko) + fiH 2 ^ H (k, ko) + 23? flWi 1 ,4(k, k ( 



(3.16) 



here 



_ 4(±)(kko) = 1± sin(2y) ±w sin(2^) +w (fW^) + ' 



and 



.A(k, k ) = 77077, 



"' 2fci2 ' ' ,u ''" V 4|k-k |i? 
'sin(2|k-ko|fl) sin(2|k+ko|i2)" 



4|k+k | J R 



4|k-ko|i2 4|k+k |i? / 

This result has a specific features therefore we put in the separate paragraph. 



Differential cross sections of electronic-vibrational transitions: now we consider differential 
cross sections of electron-vibrational transitions. Using the analytical form (see Eqs. Q2.9|) and 
( |2.10| )) of electron transition amplitude / n (k, ko,R) we can perform exact averaging over all 
directions of R. In order to make it we begin with electron transition amplitude transforming 
into infinite sum of the spherical harmonics with unit vector R as the argument 



/ n (k,ko,R) = -47r]r £ ^>T(R) [(^ +) -^ } ) jK|k + k | J R)^ m (k + k ) + 



l m=—l 



(3.17) 



(fiW + fiH) j,(|k-ko|i2)lT(k-ko) 



where the value / is even in case of r] r] n = +1 and odd if r] r] n = — 1, ji(z) = J^Ji + i/ 2 (z), and 

J 1+1/2(2) are the Bessel functions with indices I + 1/2. Plugging the transformed amplitude 
( |3.17| ) into the expression for the differential cross section (see Eq. ( J3.13Q ) and averaging over 
angular variables R using the addition theorem for spherical harmonics yields 



m=—l 



21 + 1 
An 



and £ >T(x)*7(y) 



m=—l 



21 + 1 

47T 



fl((x,y)). 



(3.18) 



Here the functions Pi(z) are Legendre polynomials [jnj. Finally we obtain the expression for 
differential cross sections of the electron-vibrational excitations 



j 



where we use following notation 



(0 
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+ 2K(^ ) ^ ) )Pi(cos^) 



(3.19) 



(0 



(|k±k |i?) TQ^) X nv ,(R)X 0v (R)dR, 



cosfi 



[k+k , k k c 



|k+k | |k-k | 

We omit the indices n,v',v in the notation g± . To some extent this expression is the gen- 
eralization of a Drukarev and Yurova formula (see Ref. M) for pure vibrational- rotational 



excitations at fixed R. However, as distinct from above mentioned result the formula (|3.19| ) 
arise from averaging over all initial and summing over final rotational molecular states so as 
only electron-vibrational excitations are taken into account. 
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Integral cross section of a pure electron transition: There exist two alternative ways to 
calculate an integral cross section (ICS). In the first one we integrate the averaged differential 
cross section (see Eq. ( p. 13 )) over all directions of the outgoing electron momentum k. Other- 
wise we can integrate the differential cross section with some fixed molecular orientation R over 
k and average over all direction of incoming electron momentum ko. Hereinafter we adhere to 
the second approach because it is more simple. The final result has the form 



crin 



sin(2kR)\ ( sm(2k RY 



2k R 



+ 



n<r> 



1 - V0Vn 



sin(2kR)\ ( sin(2fc 0J R)~ 



2kR I 



2k R 



\X Qv (R)\ 2 dR. (3.20) 



Integral cross sections for electron-vibration transitions: In order to calculate the ICSs 
for electron-vibration transitions we also adhere to the second approach, because it allows to 
achieve one's purpose without Clebsch-Gordan coefficients utilization. It's convenient before 
integration over k, k to represent the electron transition amplitudes / n (k, k ,R) (see Eqs. 
( |2.9| ), ( p.lOp ) as infinite sum over spherical harmonics with unit vectors k, and ko as arguments 



/ n (k,ko,R) =32tt 2 



-^ +) E E E E 

even I m=—l L M=—L 
k) Lq 



dl+L) 



3l (k R)j L (kR) y-(R)^ m (ko)^(R)F L M (k) + 



;(lo+L ) 



j l0 (k R) j Lo (kR) Y™° (R)C (k )r^ (R)^ (k) 



^ E E E E '■ 

odd Iq ni|) = - Iq Lq Mq = — Lq 

where summations are performed over even L and odd L in case of f]or]i = +1 and over odd L 
and even Lq in case of rjorji — — 1. In accordance with foregoing the integral cross section can 
be obtained by the following integration 

1 kn 



ainv' <— Of) = — — - dk dk 
An k 



X nv >(R) / n (k, ko, K)X 0v (R) dR 



Using the spherical harmonic orthogonality and Eq. fl3.18| ) we obtain for integral cross section 



E E(2' + i)(2£ + i) 

even I L 



a[nv <- 
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Ov) = 167T— • 

k 
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odd Iq Lq 



SloLo) 



(3.21) 



here we use the following notation 

n^ 3l (k R) 3L (k n R)X nv ,(R)X 0v (R) dR. 



(IL) 

q± 
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Figure 2: DCS for the pure X 1 S+ — > o 3 S+ transition H2 by electron impact at (a) 15 
eV and (b) 18 eV. Solid lines, presents our calculation for b = 1.35 (lower line), 1.40 
(midline), and 1.45 (upper line); full circles, SMC results of Lima, et al. [|I7J; crosses, 
MCF results of Lee, et al. ||] 

4 Applications and discussion 

Molecule H 2 , — > a 3 S+ transition: selecting the ZRPs parameters, we proceeded 
from the assumption that elements of the scattering length matrix 



(A 



a c 
c b 



are independent of number of channels N. Therefore, the parameter a can be chosen equal 
1/0.35 (see Ref. ]7j) and c, b are adjusting parameters. In our calculation, we also used the 
vibration quantum ujq = 2 x 10 -2 , anharmonicity constant aso = 5.74 x 10~ 4 and equilibrium 
internuclear distance 2Rq = 1.401. 

In Figs. H, we present our calculated differential cross sections (DCSs) for pure electron- 
impact electronic excitation for c = 0.63 and b = 1.35, 1.40, and 1.45. We also compare our 
DCSs at some selected energies with the Schwinger multichannel (SMC) results of Lima, et al. 
17] and method of continued fractions (MCF) results of Lee, et al. ]TB[. In general, there is 



good qualitative agreement. However, our and their results differs in the forward and backward 
directions, particularly for impact energies above 17-18 eV. 

Fig. |3| show our integral cross sections for pure electron-impact electronic excitation. In 
our calculation, the parameters c = 0.63, and b = 1.4 were used. We compare our ICSs with 



the SMC results of Lima, et al. [|T7J and MCF results of Lee, et al. Ps| . Our calculated ICSs 
are in good agreement, both qualitatively and quantitativly, with their theoretical data in the 
13-17 eV range. Our ICSs smaller then their results for impact energies above 17-18 eV. 
We think that possibilities of the ZRP methods allow to improve the results if 

1) to use better approximations for a potential well in which the molecule oscillates, 

2) to incorporate a separable potential into the theory (e.g., 0, |16|), 

3) account other channels of electron excitations. 

All these developments of the model could diminish the difference between the plots obtained 
by this approach, simulations "ab initio" and experiments. 



8 



12 14 16 18 20 22 24 

Impact Energy (eV) 



Figure 3: ICS for the pure — > o 3 S+ electron transition by electron impact in 
the 12-24 eV range. Solid line, present our calculation; full circles, SMC results of 
Lima, et al. []T7j] ; crosses, MCF results of Lee, et al. An experimental point is also 
shown. [[TjJ 

5 Conclusion 

We claim that the next step of our work is to unify both multichannel and higher partial modes 
descriptions. This way we evaluate amplitudes of the N2 electron-molecular scattering at low 
energies to fit recent experiments ||19|| . The results will be published elsewhere. The approach 
has one more natural generalization for multicenter scattering || |2(| 
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A The vibrational harmonics X nv (R) 



In this work we approximate the theoretical energy "curves" of molecular states by the Morse 
potentials 



U n (R) 



4a3 r , 



1 - exp -2J2/j,se n (R - R n ] 



+ U n , 



where /1 is reduced mass of molecule, u n , ae n - corresponding vibration quantum and anhar- 
monicity constant of the electronic state n, U n are energies at equilibrium internuclear distances 
2R n , and 2R - internuclear distance. The choice arbitrariness of the parameters U n can be re- 
stricted by ground state energy fixation. For that we demand ground state energy be equal 
zero. In this case 

TT ^0 . aeo 
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The vibration harmonic X nv (R) satisfied the radial Schrodinger equation 

q X nv ~\~ U n (R)X nv E nv X nv: 

o/i 

here argument R is omitted. The physically reasonable solution of this equation can be ex- 
pressed via orthogonal Laguerre polynomials L^{z) 

X nv (R) = {C nv y l l 2 z^ exp(-z/2)Li(z), 

where we use the following notation 



z = — exp(-2y'2/xa3 n (i? - R n )), 



e = ^-2r,-i. 

33™ 



The normalization factors C nv are defined by the normality condition. If functions X nv (R) are 
normalized to unity then the high accurate approximations for normalization factors are given 

T^ + v + 1) 

^nv y • 

v\^8mee n 

The energy levels E nv constitute the finite sequence 

E nv = u n (v + ^ - ae n (v + ^ + U n , v = 0,l, ..,v m . 
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